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The band structures of the periodic nanotube junctions are investigated by the effective mass 
theory and the tight binding model. The periodic junctions are constructed by introducing pairs 
of a pentagonal defect and a heptagonal defect periodically in the carbon nanotube. We treat the 
periodic junctions whose unit cell is composed by two kinds of metallic nanotubes with almost same 
radii, the ratio of which is between 0.7 and 1 . The discussed energy region is near the undoped 
Fermi level where the channel number is kept to two, so there are two bands. The energy bands are 
expressed with closed analytical forms by the effective mass theory with some assumptions, and they 
coincide well with the numerical results by the tight binding model. Differences between the two 
methods are also discussed. Origin of correspondence between the band structures and the phason 
pattern discussed in Phys. Rev. B 53, 2114, is clarified. The width of the gap and the band are in 
inverse proportion to the length of the unit cell, which is the sum of the lengths measured along the 
tube axis in each tube part and along 'radial' direction in the junction part. The degeneracy and 
repulsion between the two bands are determined only from symmetries. 
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I. INTRODUCTION 

Carbon nanotubes are one-dimensional structures 
formed by rolling up the honeycomb lattice of the mono- 
layer graphitcu Their radius and length are of nanometer 
and micrometer sizes, respectively. Qne of their inter- 
esting featureSptheoretically predictedH and investigated 
experimentally^ is that they become metallic or semi- 
conducting according to the radius and the helicity of the 
honeycomb lattice forming the tubes. Especially metallic 
nanotubes are expected to be used as electric leads with 
nanometer size. Thus we concentrate our discussion to 
the metallic nanotubes in this paper. 

A junction connecting different nanotubes can be 
formed without dangling bonds by introducing pair 
of a pentagonal defect and a heptagonal defect DI3 Such 
defects are called disclinations and are necessitated to 
form various structures composed of curved-,surface of 
graphitic layer ,BII3_for example, fullexenes,tll minimal- 
surface structurestij, torus structures^, cap stru ctures 
at the end of the tubeE-3 and helical nanotubesEalia. The 
tight binding calculations show that the electronic states 
near the undoped Fermi level are drastically changed ac- 
cording to whether the Kekule patterns started ffpp the 
different disclinations match with each other.BEJ The 
boundary between the mismatched Kekule patterns is 
called 'phason line' in Refill and the same word is used 
also in this paper. The Kekule pattern represents the 
periodicity of the Bloch states at the K and K' corner 
points in the 2D Brillouin zone of the monolayer graphite, 
at which the Fermi level exists. It means that certain 
characters of the Bloch states remain, though the Bloch 
states themselves can not be the eigen states of the sys- 
tems. In fact the electronic states of the junctions near 
the Fermi level can be described fairly well by the Bloch 
states multiplied by envelop functions as shown in this 
paper. 

The way of tiling the honeycomb pattern on the nan- 
otubes determines whether they are metallic or semicon- 
ducting. Adding only one atomic raw to a nanotube can 
change a metallic nanotube into a semiconducting one or 
vice versa. But once one knows whether the nanotube 
is semiconducting or metallic, necessary information to 
determine the electronic structures near the Fermi level 
is only about the size and the shape, i.e., one can forget 
detailed information about the honeycomb lattice such as 
the direction of the honeycomb row with respect to the 
tube axis. For example, the gap of the semiconducting 
tube is in inverse proportion to its radius because of the 
isotropic linear dispersion relation at the K and K' cor- 
ner points. Another example is the conductances of the 
junction connecting two metallic nanotubes. It is deter- 
mined almost only by the ratio of the circumferences of 
the nanotubes and the ratio \E/E C \, where E c is threi 
old energy above which the channel number increases 
Both examples are independent of the helicity of the hon- 
eycomb lattice and suggest that some continuum theory 



ignoring the atomic details but including only the size 
and shape of the systems is effective. The continuum the- 
ory describing the envelop factors of the wave functions 
is known as the effective mass theory or the k ■ p approx- 
imation. Purpose of this paper is to explain the origin of 
the correspondence between the band structures and the 
phason lines_of.the periodic nanotube junctions by using 



this theory.lixH Especially, whether the bands are degen- 



erate or avoid each other is related to the symmetries as 
discussed in section IV . Furthermore dependence of the 
band structures on the size and the shape of the system 
is discussed in detail based on this theory. 



II. EFFECTIVE MASS THEORY AND ITS 
APPLICATION TO THE SINGLE NANOTUBE 
JUNCTION 

The single junction is discussed in this section, which 
provides a basis for the discussions of the periodic junc- 
tions. The [expression is changed from that of our pre- 
vious papenla and other referencest!E!J to facilitate dis- 
cussions. First of all, we explain the Bloch state of a 
monolayer graphite forming the single wall nanotube by 
the tight binding model and relate it to the effective 
mass theory.E3 Fig. [j] shows the development map of 
the nanotube. The vector R represents the circumfer- 
ence of the tube. Two parallel lines perpendicular to R 
and parallel to the tube axis are connected with each 
other to form the tube. Here we use two pairs of the 
vectors {ei,e2} and {e x ,e y }, where e x — (e*i + e^) /V3 
and e y = e-i — e\, to represent components of vectors 
on the development map. For example, the components 
of R in Fig. |l| are represented as (RijRz) = (2,5) and 
(R x ,Ry) = (7^3/2,1.5). In this paper, we concentrate 
our discussion to the metallic nanotube, so that only the 
tube of which- R\ — R2 is an integer multiple of three 
is considered.El The four vectors have the same length 
which is about 0.25 nm and denoted by a hereafter. The 
amplitudes of the wave function at q — (qi , q2 ) and at 
q + f = (31 + 1,32 + 5) are represented by iPa(q) and 
■05 (g), respectively, with integer components q\ and (72 ■ 
The wave function -0 can be represented by the Bloch 
state as ijji(q) = e i( - kiqi+k2q ^ a tpi(0) (i = A,B). When 
the metallic nanotube is not doped, i.e. ,the 7r band 
is half filled, the Fermi energy locates at the K and 
K' corner points in the 2D Brillouin zone: the cor- 
responding wave numbers (kia,It2a) are (27r/3, — 2tt/3) 
and (— 27r/3, 27r/3), respectively. The corresponding en- 
ergy position, i.e. , the Fermi level of the undoped 
system, is taken to be zero hereafter. Near the un- 
doped Fermi level, i.e., when the wavenumber k is near 
the corner point K, the wavenumber k' measured from 
the K point, (k[a, k' 2 a) — (k\a — lixji^k^a + 27r/3), 
is small so that the phase factors can be linearized as 
exp(i/cia) = wexp(i/c' 1 a) ~ w(l + ik[a) and exp(iA)2a) = 
u> -1 exp(z/c 2 a) ~ w (1 + ik' 2 a) , where w = exp(z2-7r/3). 
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Then Schrodinger equation of a simple tight binding 
model for the Bloch state becomes 



and 



V3 



ja(k' - ik' x )ip A (q) 



(1) 



(2) 



where k' x = (k[ + k' 2 )/V3 and k y = k' 2 — k[. Here 
7(~ — 2.7eV) is the hopping integral between the nearest 
neighboring sites. In this tight binding model, only the it 
orbital is considered and mixing between the a and the 
7r orbital caused by the finite curvature is neglected. The 
solution of these equations shows the linear dispersion 
relation, 



± 



2E 
\/3~fa 



(3) 



For the one-dimensional band which intersects the K 
point, the periodic boundary condition around the cir- 
cumference is Rik[ +^2^2 = 0- From this condition, one 
can show that phase difference between A sublattice and 
B sublattice is represented by 



^B(q)/ipA(q) = ±exp(irj) 



(4) 



where 77 is the angle of R with respect to e x measured 
anti-clockwise as shown Fig. [j] . 

In the effective mass theory, the wave function is rep- 
resented by 

i/iiffl = Ff 1 '($)«;(«-**) + F? ' %w^- q ^ (% = A,B). 

(5) 



where F 



K,K' 



A,B 



W 



(91-92) w (q2-qi) 



are the envelop wave 
functions, the Bloch state wave function at the K point 
and that at the K' point, respectively. This definition 
of F's is differentnfcom that of our previous paperc3 and 
other referencesE3E3 by certain factors. The reason why 
this definition is used is that the representation of the 
time reversal operation / becomes simpler as 



I(F«,F«,F 



K' 
A : 



F, 



"'" = ((FTr,(F«'y,(F«y,(F«y) 

(6) 



so that the corresponding scattering matrix becomes 
symmetric, which is important in the discussions in sec- 
tion 



[V 



When the Fermi level Ep is close to zero, spatial vari- 
ance of the envelop function is slow compared to the lat- 
tice constant, a, so that it is a good approximation to 
take only the first order term in the Taylor expansion of 
the envelop function. This approximation corresponds to 



the replacement k' x — > —id x and kl 



id y in eq. (m) 



and eq. 



from which one obtains 



-id v 



{-idy - d x )F A 



kFf 



(7) 



(8) 



Here k = \k'\ and the isotropic linear dispersion relation 
(||) is used. The equations of the envelop wave functions 
Fab f° r t ne K' corner point can be easily obtained in a 
similar way as 



(idy + d a )Fg' = kF^ 



?K' 



(idy 



d x )Ff 



kF, 



K' 



(9) 



(10) 



Hereafter the envelop wave functions F's are often simply 
called the wave functions. 

It can be seen from eq. (||) that the Bloch state wave 
function for the one dimensional band intersecting the 
K point is represented by the envelop wave functions 



i>K± = (e-' i " /2 ,±e l " /2 ,0,0)e ±j(fc "'- , 7 ) 



(11) 



where the upper sign and the lower sign represent the 
direction of the propagating waves. Here, positive direc- 
tion is taken to be from the thicker tube to the thinner 
tube, as is shown in Fig. |^. From the propagating waves 
near the K point, the other propagating waves ipK'± are 
obtained by the time reversal operation (||) as 



K'± 



(0,0,e l " /2 ,Te" 



(12) 



Note that the direction of the propagation is reversed by 
the time reversal operation /. 

When the direction of x' = x^; (j = 5, 7) is taken to 
be parallel to the circumference of each tube as shown 
Fig. ||, k' x , is quantized as k' x ,(n) — 2im/R and k' y , is 

given by k' y ,(n) = ^/k 2 — k x '(n) 2 . Here n is an integer 
representing a number of nodes around the circumference 
and R is the circumference of the tube. When k' yl (n) is a 
real number, the channel n is open and the corresponding 
wave function is extended, otherwise the channel is closed 
and the wave function shows exponential grow or decay. 
The number of the open channel is called the channel 
number. When the Fermi energy is zero, only the channel 
n — is open, and therefore the channel number is kept 
to two irrespective of R. 

The electronic states at the Fermi energy (Ep = 0) 
govern the electron transport for the undoped system, so 
discussion in this paper is concentrated to the energy re- 
gion near zero where the channel number is kept to two. 
In order to discuss the wave function in the junction part, 
the polar coordinate (r, 9) is useful. Its relation to the 
coordinate (x,y) is the usual one, i.e. ,r = y/ 'x 2 + y 2 , 
tan0 = y/x. Fig. [2] is the development map of the 
nanotube junction where the origin of the coordinate is 
defined.Q A heptagonal defect and a pentagonal defect 
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are introduced at Ps(= Q5) and Pt{= Q7), respectively. 
Thus from now on, the indices '7' and '5' are used to 
represent the thinner and the thicker tube, respectively. 
The equilateral triangles 'AOP7Q7' and 'AOP 5 Q 5 ' with 
bases '-P7Q7' and '-P5Q5' have common apex O, which 
is taken to be the origin of the coordinate (x,y) in this 
paperLJ Then the wave function satisfies the wave equa- 
tion (z 2 d 2 + zd z + &l + z 2 )F = 0, where z = kr. The 
solution is represented by Bessel functions J m and Neu- 
mann functions N m as 



Ff = e™f 2 (z) + e-"h(z), 



F= J2 e md (c m J lml (z)+d m N M (z)) 



(13) 



The boundary condition in the junction part is if){r, 6- 
7r/3) — i/j(r, 9). It is represented by 



Mqi, Q2) = V'sfai + <72, -qi - 1) 
1pB(qu 32) = 1pA(qi +92 + 1, -qi - 1) 



(14) 



where q = (91, (72) is the position of arbitrary atoms in the 
junction part. Fig. ^ shows an example where q\ = q 2l 
i.e., the points B(2i, —i — 1) and A(2i + 1) on line 

OQ are transformed into A(i,i) and B(i,i) on line OP, 
respectively, under the rotation by 7r/3 with respect to 
the origin O. The angle 6, i.e., the direction of OQ can be 
taken arbitrary. In Fig. [|, it is taken to be parallel to the 
bond, i.e., 9 = — 7r/3, only to simplify the presentation. 

From eq. ( p~4|) and eq. the boundary conditions in 
the junction part are represented by 



F% (z,6 + tt/3) = wFg(z,6) 



Fg(z,6 + %/3) = wFf(z,6) 



F«(z,6 + 7r/3) = ^FK'(z,6) 



F^'(z,9 + 7r/3) = -F^(z,e) 
w 



(15) 



(16) 



(17) 



(18) 



Similar boundary conditions are discussed by Matsumura 
and Ando.EJ But the boundary conditions, eq. ([l^), eq. 
(|l|), eq. @ and eq. (|l|) are different fr om those of 
Matsumura and Ando by certain factors due to the dif- 
ference of the definition of F A 'g . In these equations 



and hereafter, 



i27r/3 



From eq. (J15|) an d eq. (16) 



terms in eq. (|l3|) for F% and Fg are not zero only when 
m = 31+2 (I =integer ). Because the open channel n = 
has no node along the circumference, it is better fitted 
to the components with smaller \m\ in eq. (|l3| ) than to 
those with larger \m\. So we assume that one can ne- 
glect all the terms except those with I — and I = — 1 
in eq. ( |l3| ) (Assumption I). Then the wave functions can 
be written as 



(19) 



and 



F* = e 2i0 f 2 (z)-e- i6 h(z) , (20) 

where 

f m (z) = c m J m (z) + d m N m (z) (m = l,2). (21) 



From eq. (|8j) and eq. (Mj , the other two wave functions 
J and Fjf can be derived from Fj£ and Fg as 

F? =-e»h(z) + e- a9 f 1 (z) , (22) 



Ff = e 



it) 



Hz) 



-i26 



where 



h(z) = c 1 J 2 {z) + d 1 N 2 {z), 
h{z) = c 2 Ji(z) + d 2 N x {z) , 



(23) 



(24) 



by using the recursion formula of the Bessel functions and 
Neumann functions. It is easily confirmed that eq. (22) 
and eq. (|23| ) satisfy the boundary conditions eq. (17) 
and eq. (fLq). The amplitude of the open channel in the 
tube, which is denoted by a, is obtained from eq. as 




dx^ (e 1 ^ F% ±e-^F§ ) (j = 5,7) 

(25) 



for the K point. The indices + and — mean directions 
in which the electronic waves propagate. R5 and -R7 are 
the circumferences of the thicker tube and thinner tube. 
Path of integral is the straight line PjQj, the angle of 
which is denoted by r/j with respect to x axis. Equations 
for are obtained from eq. ( pq ) by replacing ± and 
K in the r. h. s with =p and K', respectively. To simplify 
the calculation, the integrations in the above equations 
are transformed as 



dx U) -» Ri 



Qi 



'3+Vj 



de 



(26) 



If variation of the wave function along the radial di- 
rections is slow near r = Rj, this replacement can 
be allowed (Assumption II). The relation between the 
amplitudes of the open channel in each tube, Sj = 



K' „,K' 



af_), and the coefficients represent- 



ing the wave functions in the junction part, c 
' (C2, d 2 , ci, di), are summarized in the followings. 



= VRjP^MA^LikR^c , 
where M is a constant matrix given by 



(27) 
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M 
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V * 








\/3 , 
2 / 



(28) 



A(?y) is a diagonal matrix, where Ai : i = 3 = e lri and 



A 22 = A 



4,4 



P(?7) is denned by eq. (p5|) as 



P(, 7 ) 





e 2 , 


0. 


° \ 


o, 


o, 


e 2 ; 


■ V 

—e l 2 


o, 


o, 








-e 


o, 


o / 



(29) 



The matrix elements of L(z) are L\\ — L33 = Ji(z), 
L12 = L 34 = N 1 (z) 1 L 2 i = £43 = J 2 {z) and L 22 = 
L44 = N 2 {z). The other matrix elements of L(z) are 
zero. From eq. (pTj), the relation between 0-7 and a§ 
is represented by three parameters, j3 = R7/R5 (ratio 
between the radii) ,(f> = 777 — 775 (angle between the tube 
axes in the development map) and z = kR§ as 



where 



a 7+ 

«7_ 



h = h>-\ 



t2, 



1 1 



«5H 
«5- 



*2 = fr- 



cos(|</>), jsin(|< 
7.sin(1j0), cos(^ 



— cos(J<?!>), — isin(!</> 
isin(^), cos(|0) 



The factors /i+ and /i_ are represented by 



h±=~(X 12 TX 21 ) 



(30) 



(31) 



(32) 



(33) 



where 



Xi.j = y/pnziJiifizWiz) 



N^z)J 3 {z)} . (34) 



The matrix in eq. ( |30| ) is called the transfer matrix 
of the junction, denoted by T s , hereafter. We assume 
here that evanescent waves can be neglected. (Assump- 
tion III) It can be easily confirmed that T s in eq. j30| ) 
satisfies the time reversal symmetry and unitarity ( See 
Appendix I). The parameter z is related to the Fermi 
energy Ep as follows. When \k\ is near zero, channel 
number is always two independent of the radius of the 
nanotubes. But as increases, the channel number 
increases firstly in the thicker tube, when \k\ exceeds 
k c = 2tt I i?5 . Owing to the linear dispersion relation 
eq. (||), z — 2nk/k c — 2irEp/E c holds where E c is the 
threshold Fermi energy corresponding to k c . The trans- 
mission rates are calculated from the transfer matrix, and 
the conductance a is obtained by Landauer's formula as 



2l\h A 
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{£tiE = i(3/4)^' 



2 X 2 - 

h3 



} + 6 



(35) 



The obtained conductance has a remarkable feature that 
it does not depend on the angle <fi. It is consistent with 
the scaling law with the two parameters Ep/E c and 
(3 = R7/R5 in Reffj. Fig. |J shows the conductances 
calculated by the tight binding model, and those by eq. 
([55|) . Agreement between the two methods is fairly good. 
When ~ 1, a is almost constant with a value near 2 in 
units of 2e 2 /h. As (3 decreases, peak structures appear at 
E slightly below E c . The former case, (3 ~ 1, is treated 
mainly in the following sections. In this case, a is almost 
the same as that of Ep — 0, which is 



Eq. 



a = 8/{(3 3 + /T 3 + 2) . (36) 
^) reproduces well the numerical results in RefJi. 



III. BAND STRUCTURES OF THE PERIODIC 
NANOTUBE JUNCTIONS 

From now on, a, which is the length of the transla- 
tion vectors ei,e 2 , and -^7, which is the hopping inte- 
gral between the nearest neighbors multiplied by , are 
taken to be units of the length and the energy, respec- 
tively. So a and -^r7 are omitted in the following expres- 
sions. Fig. |^ shows the unit cell of the periodic junctions, 
which is composed by the two equivalent nanotube junc- 
tions, where one of the junctions is rotated by tt with 
respect to the other. The unit cell in the development 
map is determined by four vectors: the circumference of 
the thicker tube i?5 = (7715,715), that of the thinner tube 
R7 = (m?, 717) , the vector connecting the two pentagons 

in the thicker tube part, £( 5 ' = (L^\ L 2 ) and that 
connecting the two heptagons in the thinner tube part 

L^ 7 ) = {L^\ L 2 ). The components of these vectors are 
referred to (e*i, e 2 ), so that all of them are integers. In our 
discussion they can be taken arbitrary integers, so far as 
both of 777.7 — ^7 and 7775 — 775 are multiples of three. The 
transmission matrix T s is determined by eq. (|30|)~ eq. 

U = 5,7), 



(M), using Rj = \Rj\ = y ,n- -r ,oj -r »ojioj 

and cos0 = (77757777 + 775777 + \{m^n-j + m,' 7 n 5 ))/(R 5 Ri). 
The transfer matrix of the unit cell T p is obtained by com- 
bining the two identical transfer matrices of the junction, 
T s . Fig. U is a schematic view of this combination. The 
scattering occurs as = T s a§ in the left junction and 
oJ^ = T s a§ in the right junction. Here superscripts L 
and R mean the the coordinate systems, origins of which 
are the left heptagon and the right heptagon , respec- 
tively. The coordinate system L is transformed to R 
with the rotation by tt and the translation L( 7 ). There- 
fore the coordinates q L and the sublattices (^4, B) of L 
are transformed to those of R as 



q' 



^ L B {q L )^ L A {q L )) 



(37) 
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As is known from eq. ( |11| ) and eq. (|12|), the Bloch state 
wave function corresponding to cJj+ and a 7 ¥ is 



(38) 



and that corresponding to and a R _ is 



iK- q n w ±( q ?- q *) n R 



(39) 



where the upper signs and the lower signs correspond to 
K and K' , respectively. From eq. (|37|),eq. (38) and eq. 

(|39|), the relation between a L 7+ and a R j- is obtained as 



= exp(iA; 7 • L( 7 )) 
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0, -w 



-h 



,,LK 



(40) 



where I7 = L 



LK 

LK' 
*5+ 



where ^5 



^ 2^ . In the same way, 
w l \ 



exp(ik' 5 ■ L( 5 



0. 



„RK' 



a 



RK 



(41) 



= r.(5) 



(5) 



Since w = 1, Zj can take ei- 
ther values of 0, or ±1 (mod 3). In order to visualize 
the integer li, the Kekule pattern with thick bonds and 
thin bonds are drawn in the development map so that 
the pentagons and the heptagons have only thin bonds 
as shown in Fig. ^| Domain boundary caused by mis- 
match of the pattern, which is called 'the phason line' 
according to Ref.Ej, occurs in the thicker tube (thinner 
tube) only when ^5 = ±1 (lj = ±1)0 The reason for this 
correspondence is that the Bloch states at the K and K' 
corner points have the same periodicity as that of the 
Kekule pattern. Note that we do not intend here that 
the Kekule pattern represents the bond alternation. In 
this paper, the Kekule pattern is used only to show the 
periodicity of the Bloch states at the K and K' corner 
points. 

The 2x2 diagonal matrices in equations ([l0|) and ( pd| ) 
are denoted by A5 and A7, respectively. Then the transfer 



matrix of the unit cell, T p , is obtained as 



Ti, T 2 * 
To, TT 



(42) 



where 



T 1 = A^ /2 ( t i 1 A 7 i 1 - £ f 2 A 7 : i i 2 )A^ 



1/2 



T 2 = A, 



1/2 



(43) 



Above equations (ph show that t T 1 = T x , and l T 2 = 



—T 2 . It means that 



(44) 



where eq. ( |74|) in Appendix is used. If x is an eigen 
vector of T p with an eigen value A, i.e., T p x — Ax, x 
is also the eigen vector of (T p + T p 1 )/2 — Re(T p ) with 
the eigen value |(A + 1/A). When A = exp(zfc^) with 
a real value of k^ p \ is the Bloch wavenumber of the 
periodic junctions. The Bloch wavenumber k^ p ' is ob- 
tained from the eigen value of the real matrix Re(T p ) as 
cos(fc (p) ) = i(A+l/A). Furthermore Re (T p ) can be block 
diagonalized as 



1, 1 
1, -1 



/ Re(Ti), Re(T 2 ) 
\Re(T 2 ), Re(Ti) 



1 



Re(Ti+T 2 ), 

0, Re(Ti - T 2 ) 



(45) 



The eigen values of Re(T p ) can be represented by those 
of Re(7\ + T 2 ) or those of Re(7\ - T 2 ). After all, the 
dimension of the matrix which has to be diagonalized 
can be made half. Since in this paper the Fermi energy 
region is considered where the channel number is kept 

to two, the dimension of Re(Ti + T 2 ) is two. The two 

(p) (p) 

energy bands of the periodic junctions, fc_p and k_, are 
obtained as 

cos(/4 p) ) = 2 ( x " + X 22 ± 



(x n - X22) 2 + 4a& - %| x 



(46) 

where Xij = Re(Ti)y and yij — He(T 2 )ij. 

We have to discuss A and T s further to obtain the dis- 
persion relation. As for the phases of A, k[ ■ LW = ELi, 
where Li is the length of LM measured along each tube 
axis direction, i.e. , Li = |iW x Ri\/\Ri\ (i = 5,7). It 
is because k' 5 ■ R5 =0 and k' 7 ■ Rj = 0, that is to say, 
k' 5 and k' 7 are parallel to each tube axis for the metal- 
lic nanotubes. For the discussion of T s , the phases p± 
of h± in eq. (33) are defined by 
and h. 



vV Tex p(^+) 

= — R/T exp(ip-). As seen from eq. ( |35| ) 
and unitarity, T and R are the transmission rate and the 
reflection rate per channel, respectively. T determines 
the Landauer's formula conductance <r, as a — 2T in 
units of 2e 2 /h. When the tubes have almost the same 
radius, 0.7i?5 < Rj < i?5, expressions about T and 
p± become simple, so the discussions in the followings 
are concentrated on this case. In this case, the trans- 
mission rate T is almost constant value irrespective of 
z = 2nE/E c . Therefore eq. ( |36| ) is used instead of eq. 
( |35| ) in the following discussions. The transmission co- 
efficient T is almost constant near unity, i.e. , T *~ 1 
and R < 1. When 1 > R 7 /R 5 > 0.7, the range of T is 
confined in 1.0 > T > 0.76. Fig. ^ shows the quantities 
that {p±(zj + i) — p±(zj)}/(Az) vs. R7/R5, where 



7r/5, Zj = jAz, j = 0, 1, 
< z < 2?r, i.e., \E\ < E c 



■,9. They represent ^§^-, for 
(z = 2nE/E c = R 5 E). When 
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R7/R5 > 0.7, the ten plots of each Ri/R§ show almost 
the same value. It means that the phases p± are almost 
proportional to E. Fig. [7] shows that the gradients of 
p+ are very close to 1 — R7/R5, while those of p-, de- 
noted by g, are almost constant between —0.1 and —0.05. 
Therefore when R^/R^ > 0.7 and \E\ < E c , approximate 
expressions can be obtained as 



pi ~ (R 5 - Rv)E , 

P 2 ~ gR 5 E (g = -0.1 ~ -0.05). 



(47) 



The following two combinations of the gradients of 
phases determine the band structures. One is f2 + = 
L5+L7+2R5 — 2Rr and the other is = £5 — L-j + 2gR$. 
The former is the length of the unit cell, which is defined 
as the sum of the four lengths, two of which are those of 
the tubes measured along each tube axis and the other 
two are those of the junctions measured along the 'radial' 
direction. 

When there is no phason line in the thinner tube, I7 = 
0, we obtain 



cos(A4 p) ) = j; coa(n+E ± |ttZ 5 ) - ^ cos(f>_£: ± |ttZ 5 ) . 

(48) 

On the other hand, when there is no phason line in the 
thicker tube, I5 = 0, the following is obtained, 



1 



R 



cos(fc 



1 OF? O 

4 p) ) = - cos(0 + £ ± -irl 7 ) cos{n_E =p -irl 7 ) 

T 3 T 3 



(49) 



Fig. || shows the comparison of the band structures cal- 
culated by equations (f48|), ( [49] ) and the tight binding 
model for case (a) where there is no phason line, i.e., 
I5 = h = 0, and for case (b) where there are the phason 
lines in only one side of the tubes, i.e., I5 = ±1, lj = 
or I5 =0,^7 = ±1. The Fermi level comes, at the high- 
est energy of the negative energy bands B3 The results 
by the two different methods agree fairly well with each 
other, especially for the two bands nearest to the Fermi 
level. The two bands are degenerate in case (a), while 
they cross near k^ = 2tt/3 in case (b). The latter shows 
that the system becomes metallic in case (b). But there 
is a significant difference between the results by the two 
methods in case (a). The gap appears at E = in the 
tight bin ding model, while it does not appear in eq. ( f48| ) 
and eq. (|49|). Because existence of this gap determines 
whether the system is metallic or semiconducting, this 
difference can not be neglected. 

To investigate the origin of this difference, we have to 
compare the transfer matrix obtained by the tight bind- 
ing model with the one obtained by eq. (|30|)^ eq. (34). 
For the transmission rate T, agreement between the two 
methods is good, as is already seen in the preceding sec- 
tion. The ratios between the matrix elements of t\ and 
ti is almost the same in the two method, as will be dis- 
cussed in the next section. Therefore the cause of this 



difference must exist in the phase factors, p±. We do not 
get complete explanations about this discrepancy yet. 
But some features can be found as below. The phase 
factors obtained by the tight binding model are written 
as p+ = x+{Rb - Ri)E + e+ and p_ = x-gRsE + e^, 
where x± and e± represent the difference from eq. (|47|). 
Fig. H shows x+ as a function of the absolute value of 
the angle, \<f>\. The number of calculated junctions is 274 
as explained in the figure caption. The values of x+ is 
almost unity in most of the cases, but it increases when 
R7/R5 approaches unity and \(f>\ is near 7r/6. Though it 
causes a change of the gradient of the dispersion, it does 
not affect existence of the gap. The range of the ratio 
X- is about < a_ < 1, and it does not affect the exis- 
tence of the gap, either. Fig. [h] shows intercepts e± as a 
function of R5 — R-? . The intercept of p+ is almost con- 
stant with the values near 0.03tt, but that of p-, denoted 
by e_, approaches zero as R$ — R7 increases. Thus we 
speculate that the nonzero intercept e_ comes from the 
effects of discreteness of the lattice and it causes the gap 
as is shown bellow. 

The second term in eq. (^8|) and eq. (^9|) can be con- 
sidered almost constant with respect to E compared to 
the first term, when fi + ^> which is valid for the 
periodic junctions with £7,-^5 3> Iffl-Rs- Therefore the 
argument in the second cosine term can be substituted 
with 2e_. Then one can see that the r. h. s. of eq. ( |4"g| ) 
and eq. (^) become larger than unity near E = 0, and 
the gap opens at k^ p ' — 0. The width of the gap W g is 
estimated as 



W g - 2axccos(T + i?cos(2e_))/0, 



(50) 



If e_ is common, W g increases as T decreases, and is in 
inverse proportion to the length of the unit cell, f2+ . 

Next, the periodic junctions in which the both kinds of 
tubes have the phason lines is discussed. The discussions 
only for Z5Z7 = 1 are necessary, since those for Z5Z7 = — 1 
can be easily obtained from it by substitution <p with 
7r/3 — cf>. When Z7Z5 = 1, the dispersion relation obtained 
by the effective mass theory is 



cos(fc^) = Tj; cos(f2 + _E) cos((5 + ) — ^ cos(fi__E) cos(<L 



1 



R 



where 



±t^V7 (5i) 



Y = sm 2 {Sl+E) sin 2 (S+) + R 2 sin 2 (Q^E) sin 2 ((5_) 



--Rsm 2 (3<f))(cos{{n + -n_)E}cos{(n + + n_)E}) + l) , (52) 



cos((5±) 



1 3 

- =F - cos(3<^) 



(53) 



Fig. [Tl] shows an example of the band structures cal- 
culated by the tight binding model and those calculated 
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by eq. (J5l]) . Agreement between the two methods is also 
satisfactory in this case. There are two types of the gap 
at E = 0. The gap opens when Y becomes negative and 
the r. h. s. in cq. ( |5l| ) becomes a complex number. It is 
denoted by type I. The gap of the other type, denoted by 
II, occurs when the r. h. s. in eq. ( fill) becomes larger 
than unity. The gap of type / appears at nonzero k^ 
while that of type II appears at T point, i.e., few = 0, 
as is seen in Fig. [ll] (a) and (b), respectively. The band 
structures with the gap of type 77 are similar to those 
for l 5 — l 7 — shown in Fig. || (a). 

By neglecting the R 2 term in Y and expanding Y with 
respect to E, the width of the gap of type I, denoted by 
Wg , can be approximated as 




6i?(l + cos 2 (2e_))/(4- 3 cos 2 ( ^<j>j) 

(54) 



Due to the gap of type I, the HOMO band and the 
LUMO band avoid each other when k^ ^ 0, n. It is 
related to the symmetry as will be discussed in the next 
section. When </> approaches zero, Wg also approaches 
zero and the band structures become similar to those for 
I5 = 0, Z7 = ±1 or l-j = 0, I5 = ±1 which is shown in Fig. 
U (b). As 4> approaches 7r/3, the gap of type II appears 
instead of type I. Its width at = tt/3, denoted by Wg 1 , 
is evaluated for — n < 2e_ < by 



W, 
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2 arccos(T - R cos(2e_ +7r/3))/fi+ (55) 



Here the correction caused by the tight binding model 
cos(2e_) is included in eq. ( p4| ) and eq. (|55|). The gap 
width is again in inverse proportion to the length of the 
unit cell, Q+. By using Taylor expansion, arccos(l — x) ~ 
%/2~e for < x <C 1, eq. (|55| ) can be approximated further 
as 



W 11 ~ 2^i?(2 + 2cos(2e_ 



7r/3))/n H 



(56) 



(|5J) and eq. | 
as a/R = V^~- 



The common factor ^/R in eq. 

be written by using eq. ( |36| ) as V-R = VI — T 
|(1 — Rj/R$). Thus relation between the gap width 
and the ratio of the circumferences R7/R5 is almost lin- 
ear. When e_ is near zero, eq. ([34]) approaches eq. 
(|56|) as <fi approaches 7r/3. It means that even when 
(f> is close to tt/3 so that type of the gap becomes 77, 
eq. ( p4[ ) can be used approximately as the gap width. 
In Fig. [l| the solid line shows W^{cj))/wf(n/3) = 
2sin(3</>/2)/ yjA — 3 cos 2 (30/2) as a function of <fi. It can 
be seen that the gap increases as <f> approaches 7r/3. The 
dotted line shows that for Z5Z7 = — 1, which is obtained 
from the solid line for Z5Z7 = 1 by reversing the horizontal 
axis. 



IV. DISCUSSIONS BASED ON SYMMETRIES 

In this section, we discuss what can be said by consid- 
ering only the symmetry of the junctions without solving 
the effective mass equations. To discuss the symmetry, 
the scattering matrix S is used. Symmetric properties of 
S and its relation to the transfer matrix are summarized 
in the Appendix. The scattering matrix S determines 
the outgoing waves for the incoming waves as 



c*5- 

«7+ 



5, t 

t, r 7 



«54 

a f- 



(57) 



where ou 



! (a+,a+ ) and dl — (a+ , aq:) , each 
component of which is defined in eq. P^)- Note 
that the order of K and K' is reversed between 
and OL-. Consider the operation Q\ defined as 
Qi(Ff ,F*,Ff,Ff) = (-Ff, Ff ,F# , -Ff ). The 
amplitude a is transformed by this operation Q\ as 



)id± = ± 



0, 1 

1, 



Ct± — ±tTiCt± 



(58) 



Since the effective mass theory equations (H) , (P|) , (P|) , (|1Q) , 
and the boundary conditions ( 15),(fu|),(|l7j),(|18]) are in- 
variant under the operation Qi, 



( r 1 r j a 1 =r j (j = 5,7) 



(59) 



(60) 



means that Tj is antisymmetric. On the other 
is symmetric owing to time reversal symmetry. 



Eq. (|5 

hand 

So Tj is diagonal. From eq. (59), eq. (|60j) and unitarity 
of S, one can write r., and t as 



and 



Re 1 



t = V5v<»«+«t) 



1, 

0, -1 



cos(/), 
isin(/), 



(J = 5,7), 



isin(/) 
cos(/) 



(61) 



(62) 



where R and T = 1 — R are the reflection rate and the 
transmission rate, while 9j and / are certain real values. 
The meaning of / is discussed in the following. 

Consider operation Q2 shown in Fig. |l3| , where the 
thinner tube part is fixed, but the thicker tube part is 
rotated by 7r/3 in the development map. Under the op- 
eration Q2, the upper development map is transformed 
into the lower development map, where the angle of the 
circumference of the thicker tube increases by tt/3 as 
775 = 775 +7r/3. Then the angle between the two tube axes 
defined by 777 — 775 decreases by tt /3 as <p' = <j) — tt/3. The 
two development maps in Fig. [l^ correspond to an identi- 
cal junction.a The difference is only how to draw the cut- 
ting line on the honeycomb plane of the junction. There- 
fore the 5* matrix of the upper development map becomes 
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the same as that of the lower one after an unitary trans- 
formation corresponding to the operation Q2- Following 
the same discussions as those for the boundary conditions 
@, ©, © and ©, (F$,FK,FZ',FK')\ B+n/3 = 
(±Ff,wFf,wF«,±;FK)\o. By using it and V > 5 = 
775 + 7r/3 in eq. (j25|) , 

Q2Ci5± = ±iaia£± (63) 

while Q 2 c*7± = Q?7±. From this symmetry, one can derive 

r 5 ((f> - tt/3) = -<7ir 5 {<p)ai 
£(0 — tt/3) = -ioi*(0) , (64) 

which leads 

/(0) = /(O) + ^- (65) 

Lastly, we consider the coordinate transformation from 
the right-handed coordinate to the left handed co- 
ordinate Q3, (x, y) — > (~x,y). It causes r/j — > 
— 77j, — * —0, and exchange between the sublattices, 
Q 3 (F% ,Ff ,Ff ,Ff ) = {FK,FK,FK',FX'), so that 

Q 3 a± = f q ^ a ± = ^4 • (66) 

The result deduced from it is 

r,j (-</>) = o- 2 rj (0)er 2 (j = 5, 7) 
t(-<f>) = £T 2 *(0)o-2 , (67) 

which indicates that /(0) = in eq. (|65|). The trans- 
fer matrix obtained from the above discussions has the 
same form as eq. © and eq. @. Here it should be 
noted that these results can not be applied strictly to 
the scattering matrix when higher order terms of k ■ p are 
included in the effective mass equations, which make the 
equations to be variant under the operation Q%. Thus 
the invariance under the operation Q\ holds within the 
linear approximation with k' . Fig. |l4] show how the Q\ 
invariance is accurate for the S matrix calculated by the 
tight binding model. It shows values of Im(t2i/tii) as 
a functions of <f>. The range of the indices of junctions 
(7715,715) - (7717,717) and that of the energy are the same 
as those of Fig. ^ and Fig. It can be seen that these 
plots coincide with tan(30/2) quite well. The other ra- 
tios between the matrix elements obtained by the tight 
binding model also coincide with those of eq. ( |6l| ) and 
eq. (||) fairly well: |f£ - 1| < 0.1, \^ + 1\ < 0.2, 
< 0.1, and |Re(fei)|||i| < 0.2 for the 274 junctions 
inFig. |lj. 

The ratios between the matrix elements of the scat- 
tering matrix is determined in this way. But as for the 
factors h±, what can be known from the above discus- 
sions is only that they are periodic even functions of (j> 
with the period tt/3. One has to solve the effective mass 
equations to get more information on them. But there is 



an important result obtained solely from the argument 
of the symmetries. As is shown in the preceding sec- 

(v) (%>) 

tions, the two bands Izf and k_ are degenerate when 
neither of the two tubes has the phason line. On the 
other hand, the two band avoid each other when <j> ^ 
( (f> 7^ 7r/3 ) and both of the two tubes have the phason 
lines Z5Z7 = 1 = —1). From eq. (f46f), these can 

be explained as below. The origin of the degeneracy in 
the former case is that both of Re(T\) and Re(T 2 ) are 
diagonal and Re(Ti)ii = Re(Ti) 22 . The origin of the 
repulsion between the bands in the latter case is that 
Re(T 2 ) has nonzero off-diagonal elements. These origins 
can be shown only from the results in this section without 
solving the effective mass equations. 

V. SUMMARY AND CONCLUSION 

The band structures of the periodic junctions com- 
posed by the two kinds of metallic nanotubes are con- 
sidered by combining the transfer matrix of the single 
junction, T s . The transfer matrix T s is obtained in an 
analytical form by the effective mass equations with as- 
sumptions I, II and III. The region of energy is \E\ < E c , 
where E c is the threshold energy above which more than 
two channels open. By combining the T s 's, the disper- 
sion relation of the periodic junctions is also obtained 
analytically. Discussions are concentrated to the case 
when the two tubes have almost the same radius, i.e., 
0.7 < -R7/.R5 < 1. In this case, the transmission rate T 
per channel is near one and almost constant with respect 
to the energy. Agreement between the band structures by 
the effective mass theory and those by the tight binding 
model is satisfactory. Correspondence between the pha- 
son lines and the band, structures near undoped Fermi 
level discussed in Ref.ta appears naturally. When there 
is no phason line (case i), the degenerate bands appear. 
When only one of the two kinds of tubes has the phason 
lines (case ii), the two bands cross with each other near 
k = 27r/3. But in case (i), there is a significant differ- 
ence between the results of the two methods. The gap 
appears at k = in the tight binding model, while it does 
not appear in the effective mass equations. The origin of 
this difference is that the values of the phase factors p- 
at E = are different between the two methods. But ex- 
planation for it is not enough yet. In these two cases (i) 
and (ii), the angle between the two tube axes, <j>, does not 
influence the band structures. But when both kinds of 
tubes have the phason lines, the band structures depend 
on <j) as follows (case hi). When Z5Z7 = 1 (Z5Z7 = —1), the 
gap width increases and the corresponding wave number 
changes from k^ ~ 2tt/3 to fc&> = 0, as (f> increases 
from to tt/3 (decreases from tt/3 to 0). The width of 
the bands and the gap near the undoped Fermi level is in 
inverse proportion to the length of the unit cell which is 
defined as the sum of the four lengths, two of which are 
those of the tubes measured along each tube axis and the 
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other two are those of the junctions measured along the 
'radial' direction, i.e., R5 — R7. It is also found that the 
ratio between the matrix elements of the transfer matrix 
T s is determined only by the symmetries. By using only 
the symmetries and without solving the effective mass 
equations, the degeneracy of the bands for case (i) and 
the gap at k ^ for case (iii) can be derived. 

In this paper, it is found that the effective mass theory 
is very useful to analyze the electronic states near the 
undoped Fermi levels. Though the discussions is limited 
to some type of the periodic nanotube junctions, similar 
phenomena have been found in other systems. One of the 
examples is another type of the periodic nanotubes com- 
posed by only one kind of the tubeJlS Another example is 
pairs of disclinations in the monolayer graphiteB These 
electronic states can be also classified according to the 
phason line pattern. We expect that similar discussions 
can be applied also to these cases. 
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APPENDIX 



In eq. fl69|), r\ and T2 are an iVi x N\ matrix and an 
Ni x N2 matrix, respectively, and they show the reflec- 
tion rate. On the other hand, the block matrices t and 
t' show the transmission rate. The conservation of the 



probability flow \xi n \ 
tary, that is to say, 



= I ^ou.* 1 2 guarantees that S is uni- 



5 f = S~ 



(70) 



In the absence of the magnetic field, the complex con- 
jugate of the wave function represented by eq. (f^) is 
also a stationary state. It is represented by x* n — Sx* out , 
which means that 



S* = S~ 



(71) 



From eq. d70|) and eq. ( |71| ) , one can know 5 is a symme- 
try matrix t S = S, which means that n — ri, — 
and t' — l t. Here attention should paid when the base 



(<) 



wave functions are unitary transformed, x'± 

eq. 



U. 



± x ± ■ 



For arbitrary unitary matrices U± , eq. (|70|) holds in the 



representation x' . In contrast to it, eq. (71) holds only 



when U® = uf*. 

When the two leads have the same channel number, i.e. 
, N± — N2, one can consider the transfer matrix T in- 
stead of the scattering matrix S. By the transfer matrix, 
the propagating wave in the lead 2, X2 = (a?2+, ^2-), is 



determined by that in the lead 1, x\ 

x 2 = Txi , 



(xi+,xi-) as 



(72) 



We consider the junction between a left lead 1 and a 
right lead 2 in the absence of magnetic field. When there 
is no magnetic field, the Hamiltonian H can be taken to 
be real H* = H . In this case, if ip is a stationary state, 
Hip = Eip, ijf is also a stationary state, Hip* = Eip*. 
The direction of the probability flow of ip* is opposite 
to that of ip. In each lead, the propagating wave is 
represented by ^f^iix^ip^ + x^_ip^*) for the lead i 
(i = 1, 2). Here the j'th left going wave ipj is taken to 

be complex conjugate of j'th right going wave ipy . They 
are normalized so that the probability flow is represented 
by la^l 2 — la^'l 2 in the each lead. The integer AT, is the 
number of the extended state in each tube and called the 
channel number. By solving the Schrodinger equation, 
the outgoing wave x ou t =' {x\-,X2+) is determined by 
the incoming wave Xm =* {x\ + ,X2-) by the scattering 
matrix S as 



where 



T 



til, t\2 
^21) ^22 



(73) 



From the time reversal symmetry, one can show in 



1' 



ti and to 



t2 by the similar discussion 



to that about S. Conservation of the probability flow is 
represented by t\t\ — 1|*2 = 1 and H\t% — t t2t\ = 0. From 
these, the inverse matrix of T is 



T = 



-%, 



% 



(74) 



The relation between S and T is represented as 



t 2 



-(iA)n 



(75) 



•&out — Sxi 



(68) 



where 



FIG. 1. Development map of the nanotube. 



S = 



n, t 

t', r 2 



(69) 
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FIG. 2. Development map of the nanotube junctions. It is 
similar to that of RefM The lines 'EP 7 \ 'P7P5' ,'PgG' are 
connected and become identical with the lines 'FQr'/QrQs' 
and 'Q 5 H', respectively. The rectangles 'EP7Q7F' and 
'P5GHQ5' form the thinner tube and the thicker tube, re- 
spectively. 'P7P5' is the rotated 'Q7Q5' by angle of 60 de- 
grees and the quadrilateral 'P7P5Q5Q7' forms a junction part 
with a shape of a part of a cone. A heptagonal defect and a 
pentagonal defect are introduced at Pt(= Q7) and Ps(= Q5), 
respectively. Lines 'EP7' and 'FQ7' are parallel to the thinner 
tube axis and lines 'P5G' and 'Q5H' are parallel to the thicker 
tube axis. The direction of the circumferences of the tubes in 
the development map is represented by their angles 775 and 777 
measured anti-clockwise with respect to the x axis denned in 
Fig. jjj. Thus the angle between the axes of the two tubes, cf>, 
is denned as <f> = 777 — 775. 

FIG. 3. An example of rr/3 rotation. B[ and A[ indicates 
the B lattice and A lattice with their positions 2iel — (i + l)el> 
and (2i + l)el — (i + Vje^, respectively. They are aligned 
along the line OQ which is parallel to the bonds. The angle 
between the line OP and the line OQ is rr/3. The points Ai 
and Bi are aligned along the line OP and their positions are 
iel + ie2. By the rotation 7r/3 with respect to the point O, 
the points B[ and A[ are transformed into the points Ai and 
Bi, respectively. 

FIG. 4. Landauer's formula Conductances calculated by 
the tight binding model and those by the effective mass the- 
ory. The former is shown by plots and the latter is shown by 
solid lines. It can be seen that agreement between them is 
good. The horizontal axis is the energy in units of the abso- 
lute value of the hopping integral, I7I = —7. The region of 
energy is \E\ < E c ~ 0.31 |^y| , where the channel number is 
kept to two. The conductances of junctions connecting the 
(2i + 4, 2i + 1) tube and the (10,10) tube for i = 0, 1, 2, 3 are 
shown. Values of R7/R5 are attached to the corresponding 
plots. 

FIG. 5. Development map of the periodic junction. The 
positions of the pentagonal defects and heptagonal defects are 
denoted by symbols 5 and 7, respectively. The circumferences 
of the thicker tube and the thinner tube are represented by 
vectors R$ and R7. The upper bold line is connected with the 
lower bold lines so that the points connected by the circum- 
ference vectors become identical. Vectors connecting the two 
pentagonal defects and the two heptagonal defects in each 
tube part are denoted by I/ 5 ) and L( 7 ', respectively. The 
four vectors R$, R7, I/ 5 ) and I/ 7 ) determine the bond net- 
work of the periodic junction treated in this paper uniquely. 
In this figure, R 5 = (2,5), R 7 = (1,4), L? 5 ) = (4,-4) and 
£(7) = (3 ; —3). In this figure, the thicker tube has the phason 
line represented by the dotted line where the Kekule pattern 
becomes incommensurate. An area between the two dotted 
lines is the unit cell of this periodic junctions 



FIG. 6. Schematic view of combination of two equivalent 
junctions to form the unit cell of the periodic junctions. 



FIG. 7. The gradients ^f, for < z < 2tt, 
i.e., \E\ < E c (z = 2tyE/E c = R 5 E). They 
are evaluated with {p±(zj+i) — p±(zj)} / (Az) , where 
Az = tt/5, Zj — jAz, j — 0,1, •••,9. The horizontal axis 
is the ratio of the circumferences, R7/R5. The gradients 
are almost same as 1 — R7/R5 , while are almost con- 
stant between —0.1 and —0.05. As a reference, the solid line 
representing 1 — R7 /R5 is shown. 

FIG. 8. The band structures of the periodic junctions for 
one of the two kinds of tubes has no phason line. The ver- 
tical axis is the energy in units of the absolute value of the 
hopping integral, |y| = -7. (a)i? 7 = (1,10), R 5 = (3,12), 

iX 7 ) = (4,-5) and L(* 5 ) = (6,-3). R7/R5 ^ 0.766 and 
<j> ~ 0.034tt. Is — h — 0, i.e., neither of the two tubes has the 
phason line, (b) R 7 = (7,1), R 5 = (6,3), iX 7 ) = (3,-2) and 
Z/5) = (9, -3). R7/R5 ~ 0.951 and (j> ~ 0.070tt. h = -1 and 
h — 0, i.e., only the thinner tube has the phason line. 



FIG. 9. Gradients 2Z±kl at z = where 

z = kRf, = 2_Ei?5/(\/3|7i) calculated by the tight binding 
model and normalized by those obtained by the effective mass 
theory, i.e., (1 — R7/R5). They are denoted by x+ an d plot- 
ted for 274 junctions. The horizontal axis is absolute values of 
the angle \<j>\ between the two tube axes in the development 
map in units of n. The range of the indices of the tubes, 
(m, n), of the calculated junctions is 1 < min(m,n) < 8 and 
\m-n\=3i, i = ~ 3 which satisfy 1 > R T /R B > 0.7. The 
error bar indicates that maximum and minimum of the x+ for 
each junction in the energy region \E/~/\ < 0.1. That is to say, 
it represents deviation from the linearity in this energy region. 
To evaluate this error bar, the scattering matrices for 9 dif- 
ferent energies equally spaced in this energy region are calcu- 
lated about each junction. Open diamond plots, closed circle 
plots, and open square plots correspond to 1 > R7/R5 > 0.95, 
0.95 > R7/R5 > 0.9 and 0.9 > R7/R5 > 0.7, respectively. 

FIG. 10. Values of the phases p± at E — divided by rr 
calculated by the tight binding model. They correspond to 
the intercepts e± divided by rr. In the effective mass theory, 
they are zero. The range of the calculated junctions is same 
as that of Fig. ^. The horizontal axis is difference between 
the circumference of thicker tube and that of the thinner tube, 
R5 — R7, which corresponds to the length of the junction part. 
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FIG. 11. The band structures of the periodic junctions for 
both of two kinds of tubes have the phason lines. The ver- 
tical axis is the energy in units of the absolute value of the 
hopping integral, | 7 | = -7. (a)R 7 = (10,1), R 5 = (7,7), 

lX r ) = (6,-5) and = (13,-6). R 7 /R 5 ~ 0.869 and 
<t> ~ 0.14tt. h = -1 and l 5 = 1. (b) R 7 = (6,6), R 5 = (7,7), 
LOT = (4,2) and Li 5 ) = (10, -3). R7/R5 ~ 0.857 and = 
?7 = —1 and h = 1. 



FIG. 12. The dependence of the gap width on the angle <j> 
scaled by its maximum value when both of the two kinds of 
tubes have the p hason lines. It is approximately represented 
by 2 sin(3(^/2)/-^/4 — 3 cos 2 (3</>/2) as discussed in the text. 
The solid line and the dotted line correspond to the cases 
Z5Z7 = 1 and I5I7 — —1, respectively. 



FIG. 13. The operation Q2 which fixes the thinner tube 
part but rotates the thicker tube part by 7r/3 in the develop- 
ment map. The upper development map is transformed into 
the lower one under the operation Q2- These two develop- 
ment map correspond to the identical junction, denoted by 
the (2,2) - (2,5) junction. 



FIG. 14. Imaginary part of £21 /in as a function of angle 
between the two tube axes 4>. The ranges of the junctions and 
energies are same as those of Fig. [| They are fitted well by 
tan(|0) represented by the solid line. 
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